Iterated and irreducible pion-photon exchange in nuclei 

N. Kaiser 

Physik Department T39, Technische Universitat Miinchen, D-85747 Garching, Germany 



Abstract 

. . . We calculate the contribution to the nuclear energy density functional which arises 

I from iterated pion-photon exchange between nucleons. In heavy nuclei, this novel charge 

■ symmetry breaking interaction leads to an additional binding of each proton by about 

CN ■ 0.2 MeV. Compared to that the analogous effect from irreducible pion-photon exchange is 

O ', negligibly small. As a possible mechanism to resolve the Nolen-Schiffer anomaly we pro- 



pose the iteration of one-photon exchange with an attractive short-range NN-inter action. 
The corresponding energy per proton reads: E[pp\ = (2a/157r^)(7r^ — 3 -|- Qln2)Appkp 
I with pp = kp/ Svr^ the proton density and App ~ 2 fm an effective (in- medium) scattering 

length. Hints for such a value of App come from phenomenological Skyrme forces and 
from the neutron matter equation of state. 
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2 ■ A classic problem in nuclear structure theory is to understand the binding energy differences 

between mirror nuclei (i.e. nuclei with the same mass number A = Z + N but with the proton 
number Z and the neutron number interchanged). If the strong nuclear force is charge 
symmetric, then these binding energy differences can be directly related to the well-understood 
^ I Coulomb interaction between the protons. However, as shown long ago by Nolen and Schiffer 
^ ' [1, 2] the experimental binding energy differences are systematically (by about 7%) larger 
than the ones calculated with a charge symmetric strong interaction and realistic nuclear 
^ . wave functions (reproducing e.g. elastic electron-scattering data). This discrepancy which 
ranges from fairly light up to the heaviest available mirror nuclei is called the Nolen-Schiffer 
anomaly. It is generally agreed that both nuclear correlations [3] and a charge-symmetry 
breaking strong interaction are important for its understanding. 

Recently, Brown et al. [4] have performed a systematic study of these binding energy 
differences the mass region A < 60 using the Skyrme-Hartree-Fock method. They have found 
that the anomaly can be resolved by either dropping the (weakly) attractive Coulomb exchange 
term in the nuclear energy density functional (see eqs.(5,7) below) or by introducing a charge- 
symmetry breaking delta-force which splits the proton-proton and neutron-neutron effective 
s-wave interactions. The strength of the adjusted charge-symmetry breaking interaction comes 
out, however, a factor 3 to 4 larger than expected from some high-precision nucleon-nucleon 
potentials (such as AV18 [5] or CDBonn [6]). Henley and Krein [7] suggested an alternative 
explanation of the Nolen-Schiffer anomaly based on the quark substructure of nucleons in 
the nuclear medium. In their Nambu-Jona-Lasinio model the down-up-quark mass difference 
induces, driven by the partial restoration of chiral symmetry, a substantial reduction of the 
neutron-proton mass difference in the nuclear medium (see Fig. 2b in ref.[7]). Clearly, any 
decrease of the neutron-proton mass difference will help to resolve the Nolen-Schiffer anomaly, 
since the calculated binding energy differences of mirror nuclei are based on the free neutron- 
proton mass difference of 1.293 MeV. For a similar approach using the quark-meson coupling 
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model, see ref . [8] . In this context it should also be noted that most relativistic (scalar- vector) 
mean-field models for nuclear structure leave out the (attractive) Coulomb exchange term by 
default and thus circumvent the problem. 

In this work we will use chiral perturbation theory to calculate the leading order long- 
range charge-symmetry breaking effects generated by the combined pion- and photon-exchange 
between nucleons. We present analytical results for the corresponding contributions to the 
nuclear energy density functional. This energy density functional actually defines a general 
starting point for (non-relativistic) nuclear structure calculations within the self-consistent 
mean-field approximation [9] . It turns out that the three-loop in-medium diagrams of iterated 
pion-photon exchange and of irreducible pion-photon exchange contribute at different orders 
in the small momentum expansion. This allows us even to draw some conclusions about 
convergence properties. 
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Fig. 1: Left: The double line symbolizes the medium insertion T{p,q) defined by eq.(4)- Next 
shown are the two-loop one-photon exchange Hartree and Fock diagrams. Their combinatoric 
factor is 1/2. 

The starting point for the construction of an explicit nuclear energy density functional is the 
density-matrix as given by a sum over the occupied energy eigenfunctions ^i,. We restrict here 
the formulation to the proton-states since only these participate in the Coulomb interaction. 
According to Negelc and Vautherin [10] the bilocal density-matrix can be expanded in relative 
and center-of-mass coordinates, a and r, as follows: 
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where j\^^{ak^ are ordinary spherical Bessel functions. The other quantities on the right hand 
side of eq.(l) are the local proton density: 
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written in terms of a local proton Fermi-momentum kp{f), and the local proton kinetic energy 
density: 

Tpin = E vvi/t (f ) ■ vvi/.(f ) . (3) 
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The Fourier transform of the expanded density-matrix in cq.(l) with respect to both coor- 
dinates a and r defines a "medium insertion" for the inhomogeneous many-fermion system 
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[11]: 
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(4) 



The double line in Fig. 1 symbolizes this medium insertion together with the assignment of 
the out- and in-going momenta p± q/2. The momentum transfer qis provided by the Fourier 

components of the inhomogeneous (proton) distributions Pp{r) and Tp{f). Going up to second 



order in spatial gradients (i.e. 
functional reads: 



deviations from homogeneity) the protonic energy density 
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Here, the functional E[pj] denotes the energy per proton. The other strength function Ft[pp) 
is related to the effective proton mass by the relation: M*{pp) = M[l + 2MFr{pp)]^^, where 
M = 938.272 MeV denotes the free proton mass. We apply the density-matrix formalism first 
to the one- photon exchange diagrams in Fig. 1. From the Hartree diagram one regains the 
well-known classical result: 
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(6) 



with a = 1/137.036 the fine structure constant. The simplification of expression for the 
Coulomb potential in the second part of eq. (6) holds for a spherically symmetric proton density 
Pp(r'). The result of the one- photon exchange Fock diagram is also well-known [12]: 
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The derivative of the corresponding energy density with respect to the proton density: 
d{ppE[pp])/dpp = —akp/i:, is sometimes also referred to as the Coulomb exchange term [4]. 





Fig. 2: Iterated pion-photon exchange Hartree and Fock diagrams, 
factor is 1/2. 



Their combinatoric 



Next, we consider pion-loop corrections to the one-photon exchange diagrams in Fig. 1. 
There are many possible diagrams, some introduce vacuum polarization, others generate (con- 
tributions to) the proton electric form factor and another class builds up the 7r7-exchange 
nucleon-nucleon potential (see eq.(16) below). The two diagrams of iterated pion-photon ex- 
change shown in Fig. 2 are however special and distinguished by the feature of carrying an 
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energy denominator equal to the difference of nucleon kinetic energies. These two diagrams 
get enhanced by the large scale factor M and therefore they are expected to be dominant. The 
left Hartree diagram in Fig. 2 is easily seen to vanish as a consequence of a zero spin-trace, 
tr{a ■ I) = 0. The evaluation of the iterated 7r7-exchange Fock diagram with two medium 
insertions (on non-neighboring nucleon propagators) leads to the loop integral: 
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which has been regularized by introducing an infinitesimal photon mass m^. Taking the limit 
^ of the analytical expression in eq.(8) one gets the following infrared-finite result for 
the loop integral: 

We warn that by setting = in the loop integrand and performing some seemingly harmless 
manipulations one can arrive at an incorrect result which misses the t:\Q\/2 term in eq.(9). 
After further reduction of the integral over the product of two Fermi spheres of radius kp one 
obtains the following contribution to the energy per proton: 



27ra;^ -|- arctan 2x 
1 + 4x2 



(10) 



with the abbreviation u = kp/rriT^. We choose the value g„N = 13.2 for the pion-nucleon 
coupling constant and m^r = 135 MeV stands for the neutral pion mass. In order to recover 
the enhancement factor M in eq.(lO) one has to remember that the pseudovector ttNN -vertex 
carries the prefactor q^^n/'^M. The contribution to the strength function Ft{pp) reads: 



Fr{p,) = [dxx\u - x)\2,u' - Aux - 2x') 



2tix + arctan 2x 
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where we have made use of the master integral eq.(Al) in ref.[ll]. An in-medium diagram 
with three medium insertions represents Pauli-blocking effects. We obtain from the iterated 
7r7-exchange Fock diagram with three medium insertions the following contributions to the 
energy per proton: 
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and to the strength function F^(pp): 
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Here, s = xy -\- \piP- 
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x'^z'^ are auxiliary functions which 



emerge in the reduction [11] of a nine-dimensional principal value integral over the product 
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of three Fermi spheres of radius kp. The infrared finiteness of the results eqs. (12,13) deserves 
some closer inspection. In the actual calculation the logarithmic term Int is accompanied 
by an additive constant ln{m.,^ / m^) . However, when integrated over z the infrared-singular 
constant ln(m^/m^) drops out since the integrand is odd under the substitution z —z. 
It is also interesting to consider the chiral limit of vanishing pion mass m.^ = 0. In that 
case all integrals occurring in cqs. (10-13) can be solved in closed form. One finds a simple 
Pp/^-behavior of the energy per proton: 

E[p^ = ^^31n^/^ - tt^ - 61n2) ~ -l.SOMeVfm^ ■ , (14) 
and a density independent strength function Ft{pp): 

While these density dependences could also be guessed through mass dimension counting, the 
numerical coefficients are highly non-trivial. 




-0.6 ' 

0.00 0.04 0.08 0.12 

Pp [fm"'] 

Fig. 3: The energy per proton E[pp\ versus the proton density Pp — kp/Sn'^. The dashed and 
full line show the result of the I'-f-exchange and the iterated n'-f- exchange Fock diagram. The 
dashed-dotted line shows the result of irreducible TTj-exchange, magnified by a factor 10. 

In Fig. 3 we show the energy per proton E[pp] as a function of the proton density Pp — 
kp/3n'^. The dashed line shows the result E[pp] = —3akp/An of the onc-photon exchange Fock 
diagram. The full line gives the result of the iterated pion-photon exchange Fock diagram 
(evaluated with a finite pion mass of m^r = 135 MeV). One notices that both energy densities 
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lead in nuclei to an attraction between the protons. Furthermore, as it should be for a 
higher order correction, iterated 7r7-exchange comes out a factor of about two smaller than 
17-exchange. The density dependence of the full line in Fig. 3 is well approximated by the fit 
function: E[ppY'^ ^ — 1.4MeVfm^ • py^. Interestingly, it does not differ much from the result 
eq.(14) valid in the strict chiral limit = 0. 
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Fig.4: The strength function Fr{pp) versus the proton density pp = kp/Sir'^. The dashed and 
full line show the result of the 1^-exchange and the iterated ir^-exchange Fock diagram. The 
dashed-dotted line shows the result of irreducible ir^-exchange, magnified by a factor 10. 

In Fig.4 we show the strength function F^-i^pp) versus the proton density Pp = kp/37T^. 
The dashed line shows the result -Fr(Pp) = Sba/SGirkp of the one-photon exchange Fock di- 
agram and the full line gives the result of the iterated pion-photon exchange Fock diagram. 
One observes a strong cancelation between both contributions, in particular at higher proton 
densities Pp > 0.06 fm^'^. Therefore one may neglect with good reason the part of the energy 
density functional proportional to Fr(Pp) in nuclear structure calculations. 

For a numerical estimate of the Coulomb energies of nuclei we choose radially symmet- 
ric proton densities of Saxon-Woods type: Pp{r) ~ [1 -|- exp((r — R)/a)]~^, with a diffuse- 
ncss parameter of a = 0.54 fm and a mean nuclear radius of i? = 1.07 fm-A^/^ (where A 
is the mass number). The proton density Pp{r) of a given nucleus is normalized to its to- 
tal proton number Z = in drr"^ pp{r ) and the Coulomb energy is then calculated via: 
-E'coui = 47r /o°° drr'^ pp{r)E[pp{r)]. In Tab. 1 we list numerical values of the Coulomb energies 
per proton E^coui/-^ for a few selected nuclei: Ca, Fe, Sn and Pb. The direct Coulomb term 
(labeled 17-Hartrcc) gives of course the dominant repulsive contribution which increases with 
the proton number Z. The Coulomb exchange term (labeled 17-Fock) is weakly attractive 
and almost independent of the proton number Z. The same features hold for the iterated 
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7r7-exchange Fock diagram which gives numerical values of Ecovx/Z about half as large as 
those of the Coulomb exchange term. One can infer from Tab. 1 that iterated 7r7-exchange 
leads to an additional binding of each proton in a heavy nucleus by about 0.2 MeV. This 
unique long-range charge-symmetry breaking interaction does therefore not help to resolve 
the Nolen-Schiffer anomaly. On the contrary, it even further enhances the discrepancies. 



Econi/Z [MeV] 


Ca Fc Sn Pb 


17-Hartree 
17-Fock 
7r7-Fock 


3.94 4.72 7.44 10.43 
-0.37 -0.37 -0.38 -0.39 
-0.18 -0.18 -0.19 -0.20 



Tab.l: Numerical values of Coulomb energies per proton for various nuclei. The units of 
Econi/Z are MeV. 




Fig. 5: Fock diagrams of irreducible pion-photon exchange. Their combinatoric factors are 1, 
1/2, 1/2, 1, 1 and 1/2, in the order shown. 

Next, we consider irreducible pion-photon exchange which is of subleading order in the 
small momentum expansion (no large scale enhancement factor M is present). The relevant 
set of Fock diagrams contributing to the energy density functional is shown in Fig. 5. The 
corresponding Hartree diagrams (with two closed nucleon lines) vanish again due to a zero 
spin-trace. By opening two nucleon lines of the diagrams in Fig. 5 one encounters the NN- 
scattering T-matrix related to irreducible 7r7-exchange between nucleons. This T-matrix has 
recently been calculated in chiral perturbation theory by van Kolck et al. [13]. We reproduce 
independently their analytical result: 
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Here, a\^2 and ti,2 denote the spin- and isospin-operators of the two nucleons and Q is the 
momentum transfer between both nucleons. Obviously, the charge-symmetry breaking T- 
matrix in eq.(16) can only contribute to elastic — > np scattering. The corresponding 
energy density depends therefore on the product of the proton and the neutron density. For 
the sake of simplicity we assume these densities to be equal and obtain then the following 
contributions from irreducible 7r7-exchange to the energy per proton: 

m = - f'"^ '"if : (1 - + . (17) 



Sn^M^ [3 Jo u^{l + 4x2) 

and to the strength function F^(pp): 

The density dependence of these results (magnified by a factor of 10) is shown by the dashed- 
dotted lines in Figs. 3,4. One observes that the contributions from irreducible 7r7-exchange 
are negligibly small. The Coulomb energies per proton Ecoui/Z are affected at the level of a 
few keV. Such tiny repulsive corrections can of course not help to resolve the Nolen-Schiffer 
anomaly. For the sake of completeness we consider also the behavior of the expressions in 
eqs. (17,18) in the chiral limit, — > 0, and find: 

In order to resolve the Nolen-Schiffer anomaly a dynamical mechanism is needed which gen- 
erates additional repulsion for the protons but leaves the neutrons unaffected. Inspired by 
eq.(14), we propose as such a possible mechanism the iteration of a short-range attractive 
NN-interaction with the (repulsive) one-photon exchange. For its explicit evaluation consider 
the three-loop diagrams in Fig. 2 with the dashed line symbolizing now a heavy scalar meson. 
The ratio gs/tris of the scalar meson's coupling constant and mass merely serves to introduce 
an effective scattering length App = {gs/ms)'^M/i.T:. After performing all occurring integrals 
we obtain from the iteration of a short-range pp-interaction with one-photon exchange the 
following contributions to the energy per proton: 

E[Pp\ = ^ (tt^ - 3 + 6 In 2) App ~ 2.053 MeVfm • App pf' , (20) 
and to the strength function Fr{pp): 

a /32, ^ 13 



^r(Pp)-^(^yln2--j^,,. (21) 

The sign-convention is chosen here such that a positive scattering length App > corresponds 
to attraction. 

In order to compensate both the 17-exchange Fock contribution and the 7r7-exchange 
Fock contribution an effective pp-scattering length of App ~ 2 fm is required. This numerical 
estimate stems from the numbers in the second and third row of Tab. 1 in combination with 
eq.(20) and the fit function E[ppY"' ~ — 1.4McVfm2 ■ p^/^. Note that this estimated value 
of the effective pp-scattering length App ~ 2fm is an order of magnitude smaller than the 
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free ^5*0 NN-scattering length a(^So) ~ 19 fm [14] and therefore it subsumes lots of many- 
body dynamics being active at finite density. Most interestingly, phenomenological Skyrme 
forces give a strong indication for such a value of the effective pp-scattering length. As a 
matter of fact, the Skyrme force parameters ^-nd Xq [4, 9] encode an effective pp-scattering 
length via the relation App = Mto{xQ — l)/47r. Inserting the values in Tab. 1 of rcf. [15] one 
finds App = 1.96 fm, 2.40 fm and 1.83 fm for the forces SKXce, SKXm and SKX, respectively. 
Averaging over the variants MSkl-6 of ref.[16] one gets for comparison App — 1.66 fm. A 
further hint on App comes from the equation of state of pure neutron matter. The result of 
the sophisticated many-body calculation of the Urbana group [17] is well approximated for 
neutron densities pn < 0.4 fm^^ by a fourth order polynomial in the neutron Fermi momentum 
kn, namely: = 3kl/10M - ankf^/M^ + pnk^/M^, with adjusted coefficients = 1-247 

and Pn — 2.168. The coefficient of the term linear in the neutron density p„ = k^/Sn'^ can 
be translated into an effective nn-scattering length of Ann — SnUn/M ~ 2.47 fm and isospin 
symmetry implies App — Ann — 2.47 fm. Admittedly, these arguments for App ~ 2fm leave 
some loose ends but there seem to be several interesting cross relations. 

In summary, we have calculated in this work the contributions of iterated and irreducible 
pion-photon exchange to the nuclear energy density functional. These novel terms repre- 
sent unique long-range charge-symmetry breaking interactions. In heavy nuclei, iterated 
7r7-exchange leads to an additional binding of each proton by about 0.2 McV. As a conse- 
quence of that the discrepancies between calculated and experimental binding energy dif- 
ferences of mirror nuclei will be further enhanced. The subleading effects from irreducible 
7r7-exchange turn out to be negligibly small. As a possible mechanism to resolve the Nolen- 
Schiffer anomaly we propose the iteration of a short-range attractive NN-interaction with 
one-photon exchange. The corresponding energy per proton has the simple analytical form: 
E[pp] = (2Q;/157r^)(7r^ — 3 + 6 ln2) App k^. There are various hints from phenomenological 
Skyrme forces and from the neutron matter equation of state for an effective pp-scattering 
length with a value around App ~ 2 fm. In a next (more quantitative) step one should perform 
nuclear structure calculations (e.g. within the Skyrme-Hartree-Fock approach) which include 
all the charge-symmetry breaking terms derived in this work. The value of the effective scat- 
tering length App should be determined in a best fit to the binding energy differences of mirror 
nuclei, and one should check the relation App = MIqIxq — l)/47r to the Skyrme force param- 
eters to and Xq. Of course, a more microscopic interpretation/derivation of it would also be 
desirable. 
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